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ABSTRACT 

We investigate the effects of tidal interactions on ttie planetary orbits and stellar spin 
rates of the WASP-18 and WASP-19 planetary systems using a forward integration 
scheme. By fitting the resulting evolutionary tracks to the observed eccentricity, semi- 
major axis and stellar rotation rate, and to the stellar age derived from isochronal 
fitting, we are able to place constraints on the stellar and planetary reduced tidal 
quality factors, Q'^ and Q^. We find that for WASP-18, log(Q;) = 8.2ltgJ^ and 

log(Qp = 7.77^1^5, implying a system age of 0.5791^;^^^; Gyr. For WASP-19 we obtain 
values of log(Qg) = 6.47lg;g5 and log(Qp) — 6.75l^;7y, suggesting a system age of 
l.GOlg'ygGyr and a remaining lifetime of 0.0067;tg.gg6i Gyr. We investigate a range 
of evolutionary histories consistent with these results and the observed parameters 
for both systems, and find that the majority imply that the stars have been spun 
up through tidal interactions as the planets spiral towards their Roche limits. We 
examine a variety of evidence for WASP-19 A's age, both for the value above and for 
a younger age consistent with gyrochronology, and conclude that the older estimate is 
more likely to be correct. This suggests that WASP-19 b might be in the final stages 
of the spiral-in process, although we are unable to rule out the possibility that it has 
a substantial remaining lifetime. 
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1 INTRODUCTION 

There are currently 538 officially confirmed extrasolar plan- 
ets known to astronomjj^ a substantial fraction of which 
belong to the 'hot Jupiters' category of Jovian-type planets 
that orbit within 0.1 AU of their host star. Improved obser- 
vational precision has, in recent years, led to the discovery 
of extreme examples of this class of planet with orbital semi- 
major axes of the order of 0.02 AU or less(see, for example, 
Sasselov| ([2003|, [Hellier et aT] ([2009|, |Hebb et aX] poTo])). 



The Jovian masses and small orbital separations of these 
planets imply that they interact tidally with their host stars. 
Tidal interactions between stars and planets were first 



studied in the context of the Solar system ( Goldreich & Soter 



19661, and have also been extensively studied as they apply 



to binary systems (Hut [1980 1, but it is only comparatively 



recently that such studies have been extended to consider 
exoplanetary systems. Tidal interactions lead to long-term 
changes in the orbital parameters of a planet, specifically the 
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eccentricity and the semi-major axis (see e.g. Mardling & 
Lln|([2002l); [Jackson, Greenberg fc Barnes] ([2 008'); 'Ba rker fc' 
Qgilvie[(|2009[|; [Jackson, Barnes fc Greenberg (2009)). The 



end result of this process is the spiral-in of the planet to- 
wards its host until it reaches the Roche limit, where it 



undergoes mass transfer through the LI point (Gu, Lin & 



Bodenheimer 2003 1, becoming disrupted in the process. Pre- 



vious studies of tidal evolution have, for the most part, fo- 
cused on the effect of tidal interactions on these two orbital 
parameters whilst neglecting the evolution of other param- 
eters that are involved. 



Dobbs-Dixon, Lin, & Mardling (20041 demonstrated 



that the effects of tidal interactions extend to the rotation 
rates of the two bodies, but until recently little work had 
been done to investigate this aspect of tidal theory as ap- 
plied to exoplanetary systems. There is, however, mounting 
interest in this aspect of exoplanetary tides, with several 
studies investigating the large scale effect of the presence of 
exoplanets on stellar rotation tPont,|2009 [Alves, Do Nasci^ 



mento & de Medeiros 2010 Lanza 20101, and some work 



investigating tidal interactions specifically (Leconte et al 



2010h. It is well known that the rotation rate of an isolated 



star declines with age through the action of magnetic brak- 
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ing ( Weber fc Davis||1967 l, which leads asymptotically to a 
power-law decay of rotation rate with the inverse square root 
of stellar age ( Skumanich|| 1 972 1. Barnes (20071 presented a 
method for determining the ages of stars using their rotation 
periods and colours, based on a Skumanich-type magnetic 
braking law and calibrated using the solar rotation period, 
as an alternative to stellar model fitting. What appears to be 
less widely recognized, or at least acknowledged, is that the 
torques that act on both star and planet as a result of tidal 
interactions can affect this natural rotational evolution, and 
may be sufficiently strong to overwhelm it, at least for short 
periods of time. 

In this paper we simulate the evolution of hot Jupiter 
systems to investigate the effects on the orbital eccentric- 
ity, orbital separation, planetary rotation rate and stellar 
rotation rate of tidal interactions between hot Jupiters and 
their host stars. In Section[2] we set out the mathematical 
basis for our models, whilst in Section|3]we describe the com- 
putational methods that were used. Section[4] and Section[5] 
contain discussion of the application of our simulations to 
the specific cases of the WASP-18 and WASP-19 systems 
respectively. It should be noted that these systems were not 
selected at random, but were chosen owing to pre-existing 
disagreement between age estimates calculated using differ- 
ent methods. These transiting hot Jupiters also have two of 
the shortest known orbital periods, and are therefore par- 
ticularly susceptible to the influence of tidal interactions. 



2 TIDAL AND WIND EVOLUTION 



Following Eggleto n, Kiseleva, fc Hut|(j l998') , 'M ardling fc Lin 
( [2002 ) and Dobbs -Dixon, Lin, fc Mardling^ (2004| tidal en- 
ergy is dissipated within a star and planet whose spin a^ces 
are aligned with the orbital axis at rates defined by the tidal 
quality factors Q's = 3Qs/2fc* and Q'p = 3Qp/2kp, where fc* 
and kp are the tidal Love numbers of the two bodies. The 
eccentricity of the orbit evolves at a rate 



The stellar and planetary moments of inertia are determined 
by their effective squared radii of gyration, Qs and ap. 

Angular momentum is carried away from the system via 
a magnetically-channeled, thermally-driven stellar wind, at 
a rate described by a standard Weber-Davis model 



-/sKSlsMin(ris,f2) 



(5) 



where the stellar moment of inertia 7s = asMsRl. The phys- 
ical scaling of the braking rate is determined by the constant 
of proportionality k, and Q is the 'saturation' rotation rate 
above which the stellar magnetic field strength is assumed 
to become independent of the stellar rotation rate. 

We use the expressions of |Dobbs-Dixon, Lin, fc| 



Mardling ( 2004 1 to describe the tidal spin evolution of the 



planet to the tidal torque: 
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and the star under the influence of both tidal and wind 
torques: 
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The total angular momentum of the system evolves as 
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The star and planet have masses Ms and Mp and radii 
Ra and Rp respectively, and rotation rates fJs and Q.p. 
The orbital frequency is deflned by Kepler's 3rd law, — 
G(Ms + Mp)/a^. 

If the system is fully aligned, then the total angular 
momentum of the orbit and the axial rotation of the star 
and the planet, perpendicular to the orbit, is 



Ltot = ^vMsx T^^^^^^ + asMsR^s^s + apMpRlQp. (4) 



Dividing the angular momentum loss rate by the orbital 
angular momentum 



Lnrh = M„Ms 



Ga(l -e2 
Mp + Ms 



(11) 



and rearranging, we obtain the expression for the evolution 
of the orbital semi-major axis: 



e 1 — 
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.(12) 



By integrating this equation we obtain an estimate of the 
time remaining to spiral-in for a planet orbiting a slowly- 
rotating star, 

5 



remain 



2Q's Ms f 
innMp \Rs 



(13) 



giving a result almost identical to the estimate of the same 
quantity derived from a slightly different formulation by 



Levrard, Winisdoerffer, fc Chabrier (20091. The quantity 



a/ Rs is derived directly from the transit duration for a tran- 
siting planet, enabling tromain to be estimated from directly- 
observed quantities for a given value of Q's- 
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Leconte et al. (20101 recently suggested that 
parametrizing the tidal evolution equations in this 
manner, with the stellar and planetary tidal quality factors 
constant in time, is not equivalent to the more traditional 
approach of considering either a constant phase lag or a 
constant time lag of the tidal bulge. They further suggested 
that such equations, if truncated to O(e^), produce both 
qualitatively and quantitatively incorrect evolutionary 
histories for systems with e > 0.2, and promote the use 
of equations derived from Hut (19811. Despite using a 



constant tidal quality factor model the tidal equations that 
we use in this work are not truncated in this fashion, and 
thus do not suffer from the problems that [Leconte et al.| 
describe. 



3 COMPUTATIONAL METHOD 

An estimate of the value of the scaling constant for magnetic 
braking, k, is required to implement (|5|. Following the stan- 
dard Weber-Davis model, the rate of change of the rotation 
rate of an isolated star owing to magnetic braking can be 
described by 

Cl, = (14) 

Integrating this under the assumption that Q2 « ^li and 
t2 » ti gives 



2 

2t2 



(15) 



which leads to the standard power law for magnetic braking, 

Poet", (16) 

with a = 0.5 ( |Skumanic"h||1972[ ). This neglect of the initial 
conditions is justified by the observed strong convergence of 
spin rates to a narrow period-colour relation in the Hyades 
( [Radick et al.|[r987[ ), and in other clusters of similar ages. 
We set ^2 = ^s, Hyades, Calculated by scaling the Prot-(J- 
K) colour relationship found for the Coma-Berenices cluster 
by |Collier Cameron et al.| |2009) and using J and K mag- 
nitudes taken from the SIMBAD on-line data archive, to 
obtain an estimate for k at the age of the Hyades. This 
value is assumed to be constant throughout the evolution of 
the system. For simplicity, in all cases the planet's rotation 
is initially assumed to be tidally locked such that Qpfi = n, 
but rip is permitted to evolve independently thereafter. 

More recent work has found that the braking law expo- 
nent, a, diverges slightly from the ideal [Skumanich value in 



some environments (see e.g. [Collier Cameroii et al.| ( [20091 )). 
If we generalise equation |14| to raise fl to the power /3, then 
we derive 



t/a 



(17) 



a universal expression for k with (3 — ^ + 1. We calibrated 
our implementation of the magnetic braking power law using 



the stars listed in table 4 of Collier Cameron et al. (20091, 



finding that an exponent of a = 0.495 ± 0.002 gives good 
agreement between the measured rotation periods and the 
evolved stellar rotation rates calculated using our method. 
We found that there was little to differentiate the evolution- 
ary tracks produced using values of a within this range, and 
therefore adopted the central value of a = 0.495. 



Starting from a given set of initial conditions 
(to, So, n^.o, rip^o, fflo), we integrate the four equations ([T]), 
([§, (0, and (|l2l using a fourth-order Runge-Kutta scheme, 
adapted from algorithms in Press et al. ( 1992 I. The unit of 



time is 1 Gyr, and the integration is allowed to run for the ap- 
proximate main-sequence lifetime, tms — 7{Ms/Mq)~^ Gyr 
for a star of mass Ms , or until the planet reaches the Roche 
limit, as defined by Egglcton (1983). To calculate the stellar 
radius of gyration the mctallicity and mass of the star are 



used to select an appropriate table from Claret ( 2004 2005 



2006 20071; at each timestep we interpolate through this 



data using the current time to derive a value for Og. This is 
not an entirely accurate method, as the series of papers by 
|Claret| provides only tables for discrete values of metallicity 
and stellar mass that in many cases do not coincide with the 
metallicity of the systems being studied. In these cases we 
therefore use the table that most closely corresponds to our 
data. 

We use the observed values of e, Qs and a, and their 
uncertainties, to evaluate the goodness-of-fit statistic at each 
timestep according to 



+ - 



a{t)f 



.(18) 



To this we add a Bayesian prior on the stellar age t to obtain 
the statistic 



(19) 



where ioba and at are the stellar age estimated from the 



star's density and effective temperature (Sozzetti et al 



|2007[ ), and its associated uncertainty. The age of the sys- 
tem is taken to be the time of the step at which C is a 
minimum. 

This forward integration method was built into two sep- 
arate computational schemes: a grid search that carried out 
the integration for each node in a four-dimensional grid in 
ao-eo-log((5s)-log((5p) parameter space, and a Markov-chain 
Monte-Carlo (MCMC) optimisation scheme to determine 
the posterior probability distributions of log(Q!,), \og(Qp), 
eo and ao, and to refine the estimates of the initial parame- 
ters returned by the grid search. 

The best-fitting evolutionary history investigated by the 
grid search is taken to be the solution that returns the ab- 
solute minimum value of the C statistic, whilst under the 
MCMC scheme the best-fitting set of parameters are taken 
to be the median values of the respective posterior probabil- 
ity distributions, with 1 — a errors derived from the values 
that delineate the central 68.3 percent of the distribution. 
We choose this approach over the absolute minimum, as the 
latter strongly depends on the precise sampling of the pa- 
rameter space by any given Markov chain. 

It was assumed that the orbit of the planet would mono- 
tonically shrink throughout the integration, thus placing a 
lower limit on the initial semi-major axis of 
The upper limit on the same parameter was set at ao,mai = 
0.1 AU, a value which encompasses 95 percent of the cur- 
rent distribution for transiting planets. The ranges of the 
tidal quality factors were set to 4.0 ^ \og{Q'p) ^ 10.0 and 
5.0 < log((5'j,) < 10.0 such that the commonly accepted 
ranges of values (10^ - lO'' for Q'„ - 10* for Q'p ([Baraffe, 
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Chabrier fc Barman|20Tol )) were encompassed together with 
an additional section of parameter space. 

The initial eccentricity distribution for the ensemble of 
exoplanets is the subject of much discussion. A value of 
eo,max = 0.2 encompasses 93 percent of the present distribu- 
tion and limits our parameter space to the region of validity 
for our tidal equations, as defined by |Leconte et al^ ( |2010[ ). 
However the migration method by which hot Jupiters find 
themselves at such short orbital periods is still uncertain. 
One particular mechanism, Kozai scattering ( |Kozai|[l962[ ), 
can pump the eccentricity to large values, and it has been 
suggested that short period planets are captured at perias- 
tron from such highly eccentric orbits. If this is the case then 
a much higher upper limit on eo could be justified. 

We must thus take into account the initial times at 
which we will be starting our integrations. |Meibom fc Math-| 



2005 ) carried out a survey of tidal circularisation in bi- 
nary stars. They found that open clusters show a character- 
istic orbital period at which binaries with the most common 
initial eccentricity circularise (they define a circularised orbit 
as one with e = 0.01). This period was found to vary with 
the age of the cluster. However it is questionable whether 
exoplanets, owing to the much lower secondary mass, will 
tend to circularise at the same o rbital period. In fact|Hansen| 
(20101) applied the formalism oflMeibom & Mathieul (|2005|) 



to exoplanetary systems, finding that circularisation periods 
were generally much shorter than for stellar binaries. We 
therefore set eo,mai = 0.8, a value encompassing 99 percent 
of the current distribution for transiting exoplanets, to al- 
low for the possibility of eccentric orbits produced by Kozai 
oscillations. 

It is important to note that stellar and planetary obliq- 
uities greater than 0.0 are not considered by the integra- 
tion methods used herein owing to the algebraic formulation 
that underlies our computational methods. We therefore as- 
sume that the orbital and spin angular momenta vectors are 
aligned, or anti-aligned. 



3.1 Markov-chain Monte-Carlo simulation 

Under our MCMC optimisation scheme, the coordinates of 
the origin for the Markov chain are set to the centre of the 
parameter space. An exception is made if this lies within 
3— ac of the coordinates of the best-fitting initial conditions, 
in which case the origin coordinates are set to lie outwith the 
3 — (Tc contour in parameter space. The burn-in phase was 
judged to be complete when the minimum value of the test 
statistic from the current integration exceeded the median of 
all previous minimum test statistic values for the first time 
( Knutson et al.|2008 |. 



4 THE WASP-18 SYSTEM 

The hot Jupiter WASP-18 b orbits the star HD10069, an 
F6 star with an effective temperature Te/ / = 6400 ± 100 K 
and a V magnitude of 9.3, situated approximately 100 pc 
fr om the Earth. A fit to stellar isochrones using the model 
of 



(2000|) gives a stellar age of QmQtolfo Gyr, 



' Girardi et al 

which can be further constrained using the observed lithium 
abundance. [Hellier et al.| ( |2009[ ) (hereafter |H09[ ) therefore 
assign the system an age of 0.5 — 1.5 Gyr, a value that the 



Table 1. System parameters for WASP-18, taken from |H09| with 
the exception of the rotation period, which is derived from the 
observed v sin I under the assumption that the inclination of the 
stellar axis is 90°, the J-K colour, which was calculated using data 
taken from the SIMBAD on-line database, and A jTriaud et al.| 

[Mot . 



Parameter 



Value 



Units 



Rs 

J-K 
vsin I 

Prot,s 

System age 
Mp 



Rorh 

A 



1.25 ± 0.13 

0.278 ±0.032 
6400 ± 100 
10.77 ±0.04 
5.64 ±0.28 
0.630j 
10.30 ±0.( 
1.106 



-1+0.950 
-0.530 



-0.072 



Msun 
Rsun 

K 

km s~^ 
days 
Gyr 



M 



J up 



-0.054 

0.02026 ± 0.00068 
0.0092 ± 0.0028 
86.0 ±2.5 
0.94145299 ± 0.00000087 
5.0i|i 



'JUP 

AU 



days 



current work attempts to improve upon. WASP-18 b was the 
first planet confirmed to have a period of less than a day, 
orbiting its host in just 0.94 days. Recent measurements of 

(2010) 



the Rossiter- McLaughlin effect by Triaud et al. (2010) have 
shown that the system is well aligned; the assumption of 
spin-orbit alignment required for our analysis is therefore 
justified. In the absence of a measured rotation period for 
WASP-18, we assume that the inclination of the stellar axis 
is 90° such that sin / = 1.0 and Vrot = vsm I. 

An observed system age of 0.630 ± 0.530 Gyr was used 
for the calculation of the C-statistic for this system, in line 
with the age found from stellar isochrones by |H09| WASP- 
18 b has a radius of the order of that of Jupiter, but is sig- 
nificantly more massive. Lacking a means of calculating 
however, we set ctp = ajupUer as a reasonable estimate. 

In light of the observed system age being consistent with 
the age of the Hyades, we compared Prot,s from Table[T]to 
Prot,s{t = t Hyades), as Calculated for our derivation of k for 
the system, finding the observed stellar rotation period to 
be shorter than the expected rotation period of 7.00 days 
at the age of the Hyades. We therefore set to — 150 Gyr~ 
iM35, calculating Qs,o by scaling the period colour relation 



of Collier Cameron et al. ( 2009 1 to this age. 

The C-statistic values returned by the grid search were 
converted into probability densities according to 

P{C) = e-^, (20) 

and plotted as a set of projection maps (Fig.[T]). These in- 
dicate regions of high probability, and therefore low points 
in the C-statistic surface, for each pair of initial parameters, 
and appear to show that there are a range of values for each 
parameter that could produce plausible evolutionary histo- 
ries with a good quality of fit. This is particularly noticeable 
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Figure 1. Projection maps of tiie test statistic probability density for WASP-18. The projections indicate distinct ranges of log(Q^), 
log(Qp) and ao that give high probability, corresponding to a low value of the C-statistic. eg provides a less clear cut high probability 
region, with several individual values giving good results. The best-fitting initial parameters were found to be ag = 0.0427 AU, eg = 0.245, 
log(Q'j,) = 7.00 and log(Qp) = 7.50, giving a system age of 0.654 Gyr. 



in eo and log(Qp)', indicating that these parameters have a 
lesser impact on the orbital and spin evolution of the system 
than ao and \og{Q'g). This is, to a certain extent, unsurpris- 
ing; the current orbital eccentricity for WASP-18 b is more 
uncertain than the orbital separation, and so will have less 
influence on the value of the C-statistic. It is also interest- 
ing to note the generally sharp transitions in Fig. [T] between 
regions with P{C) > 0.4 and those with P{C) ^ 0. This 
indicates that although the values of the initial parameters 
are somewhat uncertain, they are generally well constrained 
to a smaller region of parameter space than we have allowed 
for. 

The best-fitting combination of initial parameters in 
these projection maps was found to be ao = 0.0428 AU, 
eo = 0.245, log(Q^) = 7.00 and log(Qp) = 7.50, indicat- 
ing a system age of 0.654 Gyr in good agreement with the 
existing estimate derived from stellar model fitting. Using 
the period colour relation of [Collier Cameron et al.| 12009] ) 
we calculate an age of 0.387 ± 0.024 Gyr for WASP-18 using 
the derived rotation period, inconsistent with both the age 
from isochrone analysis and the result of our grid search. 
From the evolutionary track we derive an extremely short 



remaining lifetime of 0.006 Gyr, implying that the planet is 
currently spiralling in towards its host and is on the verge of 
reaching the Roche limit. This seems incredibly short, and 
would mean that we have managed to observe WASP-18 b in 
a very short window of opportunity. This result is, however, 
at oodds with the results that we have obtained from our 
second integration method, the MCMC algorithm. 

The posterior probability distributions produced by the 
MCMC integration scheme (Fig.[2| indicate significant cor- 
relation between many of the pairs of initial orbital parame- 
ters, in particular ao and log((3s) although it is also notice- 
able in several of the other distributions. Some correlation 
is to be expected; the orbital circularisation time-scale de- 
pends on energy dissipation within both the star and the 
planet ( |Miralda-Escude|2002[ ) , and so a correlation between 
the eccentricity and the quality factors that govern tidal en- 
ergy dissipation is unsurprising. The correlation between ao 
and eo may arise as a result of the strong coupling between 
eccentricity and separation noted by [Jackson, Greenberg fc] 
Barnes (20081. The rapidity with which the orbital separa- 



tion can decrease is intrinsically linked to the efficiency with 
which energy is dissipated within the system, a process that 



6 D. J. A. Brown et al 




6 8 '0 




a 10 4 E a 10 




Figure 2. The posterior probability distributions of each pair of jump parameters produced by the MCMC analysis of the WASP-18 
system. Note the clear correlations between several of the integration parameter pairs; these arise from the dependence of tidal energy 
dissipation on the orbital eccentricity, and from the strong coupling that exists between the orbital eccentricity and the orbital separation. 
The location of the best-fitting parameters in the parameter space explored are denoted by the dashed lines. 



is gove rned by the tidal quality factors. Moreover inspection 
of ( 12 I indicates that a oc , so the form of the correlation 



between ao and \og(Q's) that is observed is also as expected. 
In fact it is surprising that little to no correlation is observed 
between ao and \og(Q'p) 

The best-fitting parameters found using the MCMC 
code are set out in Table|2] and imply a stellar age of 
0.579125^0°^ Gyr, in agreement with the existing age obtained 
through isochrone fitting, and the age derived from our grid 
search results. Fig. [3] displays the evolution of the stellar ro- 
tation for the set of best-fitting parameters; it is important 
to recall that the age estimate is not evaluated solely on 
the basis of this parameter, and that the observed orbital 
eccentricity and semi-major axis play a part as well. 

The value of log(Qs) returned by the MCMC-derived in- 
tegration scheme does not agree with the best-fitting value 
obtained from the grid search of the defined parameter 
space. Nor do the values of oq and eo. We attribute this 
to the adoption of the set of median parameters from the 
MCMC scheme to avoid chance encounters with local min- 
ima, as well as the discrete nature of the grid search. Com- 
parison of the solution with the lowest value of the C test 



statistic shows a value of log(Qs) that fits more closely to 
the grid search result, albeit still in disagreement, although 
we note that the minimum test statistic value for the grid 
search is an order of magnitude lower than that for the 
MCMC exploration of the available parameter space. We 
also note that although Figs. [T] and [2] are broadly similar in 
form, they differ somewhat in detail. The same correlations 
between parameters are visible in both figures, but the cut- 
off in log((3s) occurs at a slightly greater value in Fig[2]than 
Fig.[l] Additionally, the MCMC algorithm explores a more 
narrow range of both ao and eo parameter space than the 
grid search, but does explore the lower end of \og(Qp) space, 
a region which the grid search suggests gives poor results. 
We again attribute this to the differences in the manner in 
which the two algorithms explore the parameter space that 
we have delineated. 

Using ( 13 1 with the data from Table[2] and the best- 
fitting value of log(Q^) from the MCMC scheme, we estimate 
that WASP-18 b will reach its Roche limit 0.086 Gyr from 
now. The evolutionary track produced using the MCMC re- 
sults implies a remaining lifetime of O.OTGIq'qm Gyr for the 
planet, consistent with this value. These are, respectively. 
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Table 2. The initial orbital parameters and tidal quality fac- 
tors of the best-fitting tidal evolution histories produced by the 
grid search and MCMC integration schemes for the WASP-18 
system. The I — a error bars for the MCMC derived values are 
estimated from the parameter values that encompass the cen- 
tral 68.3 percent of the final parameter distributions, and in some 
cases are inflated by the presence of short 'tails' in the distribu- 
tions. 



Parameter Grid search MCMC Units 



ao 


0.0427 




AU 


eo 


0.245 






log(Q'J 


7.00 


o r,-, +0.90 

«-^-'--0.52 




log(Q;) 


7.50 






age 


0.654 


n ^7n+0'305 
U.O/y_Q 250 


Gyr 


^remain 


0.006 


076+°-''3° 


Gyr 



2.40 and 2.12_j 23 percent of the expected main sequence 
lifetime of the host star, which we estimate to be 3.58 Gyr. 
These remaining lifetimes strongly imply that the WASP-18 
planetary system has an short life expectancy, as found by 
|H09[ Although still quite short, these values are far more 
reasonable than the value of 6Myr derived from the grid 
search. 

The evolution of the stellar rotation period (Fig.[3f for 
the set of initial parameters returned by the MCMC inte- 
gration scheme indicates that the rotational evolution of 
the star under the influence of tides initially differs little 
from our simulation in which the rotation period is governed 
purely by a Skumariich-type magnetic braking law. However 
the evolution of the rotation period rapidly diverges from 
this ideal case, with the rate of spin-down slowing gradually 
as the star ages and the planet migrates towards its host 
under the influence of tidal interactions. The final spiral-in 
of the planetary companion causes a rapid and substantial 
spin-up of the star; during this process the stellar rotation 
period is reduced by almost a factor of two from its maxi- 
mum 6.26 day period to 3.16 days at the Roche limit, sug- 
gesting significant spin up of the host star. Furthermore, 
Fig.[3] clearly shows that, for the age estimate returned by 
our MCMC integration scheme, the observed rotation pe- 
riod of the host star is quite inconsistent with the rotation 
period of 7.96 days that would be expected at the same age 
if tidal interactions played no part in the evolution of stellar 
rotation. 

Although the best-fitting parameters that we have 
adopted show WASP-18 b in the process of its final spiral 
in, the true picture is slightly more ambiguous, as evidenced 
by the large upper uncertainty on the remaining lifetime. 
Within the 1 — a parameter ranges there exist combinations 
of the initial parameters that produce evolutionary tracks 
of a different form to that displayed in Fig.jS] These evo- 
lutionary tracks follow the stellar rotation curve produced 
by a purely magnetic braking scenario for a much longer 
period of time, and exhibit little spin up of the host star. 
Furthermore, in semi-major axis space they indicate that the 
planet stays at approximately the same semi-major axis for 
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Figure 3. The evolution of the rotation period of WASP-18 from 
the best-fitting initial conditions found by the MCMC exploration 
of the allowed parameter space (solid line), along with the evolu- 
tion that would occur in the absence of tidal interactions (dashed 
line). For reference, the observed parameters are plotted with 
associated errors, and the best-fitting age given by the MCMC 
integration scheme is denoted (dotted line). The tide-governed 
evolution diverges from the ideal rotational evolution early on, 
with the effect of tidal interactions on the the rotation of the host 
star increasing with time until the rotation period dramatically 
decreases during the final spiral-in of the planetary body. 

the duration of the star's main sequence lifetime. In these 
cases the age of the system lies towards the lower end of the 
adopted range, owing to the necessity of conforming to the 
calculated stellar rotation period, and the planet has a long 
remaining lifetime. We note however that these solutions 
form only a minority of the evolutionary tracks consistent 
with the parameter ranges that we adopt from the MCMC 
results. More prevalent were tracks of the same form as that 
displayed in Fig.[3] but with the spin-up taking place much 
later such that the stellar rotation implies an age consistent 
with the gyrochronological estimate and magnetic-braking 
only scenario. 

We attempted to further constrain the range of possible 
tidal quality factors by visually fitting to the observed pa- 
rameters and their associated 1 — a uncertainties. Starting 
from the best-fitting initial orbital separation and eccentric- 
ity, we investigated the evolution of the system across the 
range of values for both \og{Q'g) and log(Qp). We found that 
changing the value oilog(Q'p) makes no difference to the evo- 
lution of the orbital separation or rotation period, but has 
a strong effect on the eccentricity evolution of the planetary 
orbit. log{Q'g), in contrast, strongly affects the evolution of 
all integration parameters. Moreover we found that chang- 
ing the initial semi-major axis, within the adopted range, 
made a significant difference to the evolution of the stel- 
lar rotation period, whilst modifying eo merely affected the 
evolution of the orbital eccentricity. We therefore conclude 
that the two most important parameters with respect to the 
stellar rotation are Q'^ and ao. 

We were unable to constrain the permissible ranges of 
the two tidal quality factors any further, owing to the range 
of values that ao was able take, but found that for the 
majority of the possible combinations the rotational evo- 
lution of WASP-18 A gradually diverges from that expected 



8 D. J. A. Brown et al 



a) 


III 




1 1 1 1 1 




, , 1 







0.2 


0.5 


1 2 
Time /Gyr 


5 


-b) 




1 1 1 1 








i 

i 
j 

i 
i 






I.I 




0.2 


0.5 


1 2 


5. 



Time /Gyr 

Figure 4. Evolutionary tracks for WASP-18 produced using sev- 
eral different values of log(Q^) within the range consistent with 
the observed system parameters. The tracks were calculated using 
the values of eo, ao ^'Hd log((3p) adopted from the MCMC solution 
and, from left to right (at the end of the track): log(Q^j) = 8.02, 
log(Q^) = 8.21, log(Q'J = 8.70 and log(Qi) = 9.11. C/pper •panel: 
The evolution of stellar rotation period with time for WASP-18. 
The dashed line shows the evolution expected from a Skumanich- 
type magnetic braking law. Changing the value of Q'^ has no effect 
on the evolution of the rotation period with time. Lower panel: 
The evolution of the C test statistic with time for the tracks in 
the upper panel. The time at which C is a minimum is taken to 
be the age of the system given the set of initial parameters used 
for that evolutionary track. 



of a purely magnetic braking scenario, and that the star is 
eventually spun-up by a substantial amount during the fi- 
nal in-spiral of its planetary companion. Fig.|4a| shows the 
rotational period evolution resulting from several values of 
\og{Qs) within the permissible range, with \og{Q'p), ao and 
eo set to values adopted from the MCMC search. 



5 THE WASP-19 SYSTEM 

The transiting hot Jupiter WASP-19 b orbits a late G-dwarf 
star with an orbital period of 0.7888399 ± 0.0000008 days. 
The host star has been measured to have T^ff = 5500 ± 
100 K, periodic sinusoidal fiux variations that indicate a 



Table 3. System parameters and 1 — a limits for WASP-19, taken 
from the free eccentricity fit of |H10| with the exceptions of the J-K 
colour, which is derived from data taken from SIMBAD, and the 
rotation period (Collier Cameron, priv. comm.). 



Parameter 


Value 


Units 


Ms 


0.95 ±0.10 


Msun 


Rs 


93+0 05 
'^•'^■^-0.04 


Rsun 


J-K 


0.43 






5500 ± 100 


K 


I) sin / 


4.0 ±0.2 


kms~^ 


Prot 


10.5 ±0.2 


days 


System age 


^+9.0 
^•''-4.5 


Gyr 


Mp 


1.14 ± 0.07 


^^jup 


Rp 


1.28 ±0.07 


Rjup 


a 


0164+°-'"'05 


AU 


e 
i 


80.8 ±0.8 




Rorh 


0.7888399 ± 0.0000008 


days 



detectable level of intrinsic variability and activity, and a 
metallicity slightly greater than solar at [M/H] = 0.1 ± 

H10[ ). Measurement 



0.1 dex ( |Hebb et alT]|2010[ ) (hereafter 

of the Rossiter-McLaughlin angle for the system indicates 
that it is aligned ( Hellier et al.|20TT l. 

The age of the system is still somewhat uncertain, with 
|H10| able to determine a constraint of > 1.0 Gyr. We note 



that several studies (e.g. Hansen (20101; Weidner & Horne 
(| 2010[ )), as well as the exoplanet encyclopedia, cite an age 
of 0.6 ± 0.5 Gyr for the star and attribute it to |H10[ W e be- 
lieve that this originates from the abstract of |H10| which 
quotes that age as one of two possibilities. However the 
text of that paper favours the older age constraint. For 
our analysis we used an age of tage = 5.5 ± 4.5 Gyr, con- 
sistent with the isochronal fit and lower bound on the age 
quoted by |H10| and set Up = ajupiter = 0.26401. Compar- 
ing the stellar rotation period from Table[3]to the rotation 
period expected at the age of the Hyades, we found that 

Prot,s > Prot, Hyades = 8.60 dayS and thus set to = tuyades 
and Q,s,0 ~ ^s, Hyades- 

Fig. [5] shows the C-statistic probability density projec- 
tion maps produced from the grid search results for the 
WASP-19 system. They indicate a fairly small region of high 
probability for log{Q's) and log(Qp), but appear to show that 
there are broader ranges of both eo and ao that give strong 
solutions. The absolute maximum probability, and thus the 
minimum value of the C statistic, was found to occur when 
the initial conditions were ao = 0.0939 AU, eo = 0.735, 
log(Qs) = 6.25 and log{Qp) = 10.0, giving a system age 
of 2.11 Gyr that is consistent with the existing constraint of 
> 1.0 Gyr from|H10[ Using the period-colour relation of |Col-| 



lier Cameron et al. ( 2009 1 we calculate a gyrochronological 
system age of 0.899±0.037Gyr, immediately suggesting that 
gyrochronological analysis is inappropriate for this system. 

Comparing the sets of projection maps for the two sys- 
tems, it is apparent that there are similarities in the forms 
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Figure 5. C-statistic probability density projection maps for WASP-19b. A small region of high probability is clearly visible in log(Q(,) 
and log{Qp) parameter space, but in ag and eo the range of values that give high probability solutions are much broader. The point 
of maximum probability, indicating the minimum in the C-statistic surface, was found to occur at parameter space coordinates of 
ao = 0.0939 AU, eo = 0.735, log(Q'J = 6.25 and log(Qp = 10.0, giving a system age of 2.11 Gyr. 



of Figs.jTJandjs] particularly in the ao — log{Q'g) and eo — ao 
maps. However the range of probability values covered by 
the greyscale is slightly narrower for WASP-19, where the 
maximum probability is approximately 0.8, compared to 
WASP-18, where the maximum probability is almost 1.0. 
From this we deduce that the best-fitting orbital solution for 
WASP-19 is less certain than our best-fitting grid search so- 
lution for WASP-18. This lower maximum probability, cou- 
pled with the form of the projection maps, also suggests 
that the range of tidal quality factors for which valid or- 
bital solutions exist is greater for this system than it was for 
WASP-18. 

The posterior probability distributions produced by the 
MCMC algorithm exhibit a strong correlation between ao 
and log(Qs), of a similar form to that observed for the 
WASP-18 system albeit extended to lower values of log{Q'g). 
Also present is a correlation between eo and ao that is some- 
what similar to that apparent in Fig.[2j but the other corre- 
lations that were present in the distributions for WASP-18 
are absent for WASP-19. The region of parameter space ex- 
plored by the MCMC algorithm is also substantially greater 
for the WASP-19 system than it was for the WASP-18 sys- 



tem, although this is not unexpected after the comparative 
appearance of Figs. [T] and |5] 

It is also interesting to note the relative lack of agree- 
ment between the results from the grid search and the 
MCMC search methods for this system; the only search pa- 
rameter for which the two sets of results agree is log((3!,). 
This is also apparent from inspection of Figs[5]and|6] as the 
region of greatest probability in the former does not match 
the location of the adopted parameters in the latter. As for 
the WASP-18 system, we attribute this to the adoption of 
the set of median parameters from the MCMC scheme; com- 
parison of the solution with the lowest value of the C test 
statistic shows that it more closely fits with the location of 
the high probability region in Fig. [5] Fortunately the stel- 
lar tidal quality factor, the most important parameter when 
considering the evolution of the stellar rotation, is relatively 
well constrained in our MCMC-derived solution. 

The probability distribution for the initial orbital eccen- 
tricity of the WASP-19 system also shows weaker clustering 
than for the WASP-18 system, and the clustering that is 
present encompasses a much smaller range of eccentricity 
values above eo — 0. This may be symptomatic of the fact 
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Figure 6. The posterior probability distributions of each pair of integration parameters produced by the MCMC analysis of the WASP- 
19 system. The correlations between the parameters are substantially different than those for the WASP-18 system, with that between 
log((3'j,) and ag being the only particularly noticeable correlation. The location of the best-fitting parameters in the parameter space 
explored are denoted by the dashed lines. 



that the eccentricity of the system is not well known; pre- 
vious analysis of the parameters of the system found little 
difference in the fit to the observed transit lightcurve be- 
tween the cases in which the orbit was forced to be circular, 
and in which the eccentricity was allowed to fioat |H10| . In 
the latter case, the best-fitting eccentricity value was only 
0.02^0 01, and it is this value that we used in our C statistic 
calculations. 

The best-fitting parameters given by the MCMC inte- 
gration scheme, set out in Table|4] imply a stellar age of 
LGOlg yg Gyr, in broad agreement with the loose constraint 
of age > 1.0 Gyr found by |H10[ The gyrochronological age 
that we have calculated for the system also agrees with 
this age estimate but lies at the lower end of the range, 
suggesting that gyrochronology provides a possible, if un- 
likely estimate for the age of WASP-19 A. The data in Ta- 
ble[4] further imply a remaining lifetime for the planet of 
O.OOGTlJ Joel Gyr, a mere 0.4 percent of the estimated sys- 
tem age. This strongly suggests that WASP-19 b is in the 
final, spiral-in stage of its orbital evolution, a conclusion 
supported by the stellar rotation period evolutionary track 



that results from integration of the best-fitting parameters 
(Fig.[7|. As with the WASP-18 system we attempted to fur- 
ther constrain the range of possible tidal quality factors by 
visually fitting to the observed parameters; in this case we 
were able to reduce the upper limit on log(Qg) and raise the 
lower limit on qq. The updated values are given in Table|4] 
Updating the limits on these parameters have no effect on 
the 1 — a limits of the age or remaining lifetime. 

The observed spin rate of WASP-19 b is inconsistent 
with evolution governed only by magnetic braking; Fig.jT] 
shows that, at the system age implied by our MCMC solu- 
tion, the observed period of 10.5 ± 0.2 days is substantially 
less than the period of 15.9 days implied by our magnetic 
braking only model. It is possible that WASP-19 b is in the 
process of spiraling-in to the Roche limit, spinning up its 
host as it does so. 

Fig. [8] shows the rotational period evolution resulting 
from log(Qs) values within the 1 — a limits returned by the 
MCMC search scheme. These reinforce the idea that WASP- 
19 b is on the verge of reaching its Roche limit; for the best- 
fitting tidal quality factors the end of the evolutionary track 
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Table 4. The initial orbital parameters and tidal quality fac- 
tors of the best-fitting tidal evolution histories produced by the 
grid search and MCMC integration schemes for the WASP-19 
system. The I — a error bars for the MCMC derived values are 
estimated from the parameter values that encompass the central 
68.3 percent of the final parameter distributions. 



Parameter Grid search 



MCMC 



Updated 
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ao 


0.0939 


n (T^l 7+0.0228 

U.UJl(_g Qj^4g 


eo 


0.735 


0017+° °'597 


log(Q'J 


6.25 


-0.95 


log(Q;) 


10.0 


6.75li;«? 


age 


2.11 


1 60+^'*'* 


^remain 


0.0038 


0-006711;;°^? 



n n^i '7"t"0.0228 
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6 47+0-6'' 
'^•^'-0.95 



AU 
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Figure 7. The evolution of the rotation period of WASP-19 from 
the best-fitting initial conditions found by the MCMC exploration 
of the allowed parameter space, along with the evolution that 
would occur in the absence of tidal interactions (dashed line). 
The tide-governed evolution follows the gradual spin-down of the 
isolated evolution closely until spiral-in of the planet begins, at 
which point the stellar rotation period rapidly decreases. The 
observed rotation period is plotted with associated errors, and 
implies that WASP-19 b is in the final spiral-in stage of its orbital 
evolution. Furthermore it is irreconcilable with the expected spin 
down from magnetic braking alone, and can only be explained by 
invoking spin-up during tidal interactions. 



is at Prot ~ 10.5 days compared to the observed P^ot ~ 
10.5 ± 0.2 day, placing the planet precisely on the Roche 
limit and suggesting that we have been fortunate to observe 
WASP-19 b at all. 

Much as with the WASP- 18 system though, this is not 
the entire picture presented by our results. The upper limit 
on tremain implied by our adopted solution is 1.114 Gyr, a 
significant 69.6 percent of the estimated age and 13.7 percent 
of the MS lifetime of the host star. This strongly suggests 
that there exist solutions that more closely follow the stellar 
rotation evolution expected when only magnetic braking is 
acting on the star. To investigate this we characterised the 




Figure 8. The evolution of stellar rotation period with time for 
WASP-19 for a range of log(Q'j,) values. The dashed line shows 
the evolution expected from a Skumanich-type magnetic braking 
law, whilst the three solid lines are the evolution that results 
from the best-fitting eg, ao a^nd log(Qp) given by the MCMC 
algorithm, and, from left to right (at the observed system age), 
log(Q^) = 6.14, log((3'J = 6.47 and log(Q'J = 6.52. These tracks 
imply that WASP-19 b is very close to the Roche limit, and will 
reach it in less than 10 Myr. 



orbital evolutionary tracks for each \og{Q'g) — ao node from 
the grid search that was consistent with our results from the 
MCMC algorithm. Since our investigation into WASP-18 b's 
future evolution led us to the conclusion that eo and \og{Q'p) 
have little effect on the evolution of the stellar rotation, we 
set these parameters to the MCMC values from table|4] 

We found that for each value of ao there was only a 
very narrow range of log(Qs) values that gave solutions with 
an age 55 1.0 Gyr and a rotation period at that age consis- 
tent with the observed value. In fact only 13 percent of the 
nodes investigated gave such solutions. The majority of the 
nodes gave solutions that either did not reach an age of 
1.0 Gyr, instead showing the distinctive spiral-in signature 
at younger ages, or that followed the magnetic braking-only 
track for much longer such that the spin-up induced by the 
final spiral-in was insufficient to reduce the stellar rotation 
period back to the observed value. Tracks that fell into this 
latter category often returned ages consistent with our gy- 
rochronological calculation, although some did return older 
ages consistent with the 1.0 Gyr lower limit owing to the 
influence of the orbital eccentricity and semi-major axis. 

This raises questions concerning WASP-19 A's true age. 
There is evidence from isochrone fitting that the star is 
older than 1 Gyr. This is supported by the investigation 
of [HlOl into the space velocity of the star, which found 
that 65 percent of simulated stars with similar stellar prop- 
erties in a small volume around WASP-19 A were older 
than 1.0 Gyr. But the rotation period of WASP-19 A im- 
plies a younger age of 0.899 Gyr through gyrochronology. 



The lithium abundance of log(y4[Lz]) < 1.0 quoted by HIO 
while consistent with either estimate as it only constrains 
the age to greater than 0.6 Gyr, would appear to support a 
younger age, as the presence of any lithium tends to rule out 
a more evolved star. There therefore appear to be two main 
possibilities; either the star is old and has been spun up by 
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the infall of the planet, or the star is younger, still following 
its natural spin down and more dense and/or hotter than 
expected for its age. Both of these are consistent with the 
results that we have obtained from our search methods, so 
which is the more plausible? 

5.1 Is W19 young or old? 

We consider our results for WASP- 19 in the context of a 
larger ensemble of exoplanetary systems with similar stel- 
lar and planetary properties, assuming that the value of 
log(Qs) = 6.47 adopted for WASP-19 is applicable to the 
other systems under consideration. We find that for 7 of the 

9 planetary systems in the sample the remaining lifetime 
is < 1 Gyr. Furthermore we find that WASP-4 b has a re- 
maining lifetime of 0.017 Gyr, and calculate that it lies at a 
distance of only 1.4 times the Roche limit from its host star; 
this places it in apparently similar circumstances to WASP- 
19 b. WASP-2b and WASP-lOb have remaining hfetimes at 
least one order of magnitude greater than those of the rest 
of the sample, but they also have significantly longer orbital 
periods. 

It is interesting to consider the range of the stellar 
tidal quality factor that produces certain threshold values 
of tremain- We find that, if we exclude the longer period 
WASP-2 and WASP-10 systems, log(Q;) > 7.10 leads to 
an estimated tremain ^ 1.0 Gyr for all planets in the sam- 
ple, whilst a remaining lifetime greater than 0.1 Gyr requires 
log(Q^) > 6.12. For WASPs -2 and -10 we find that the val- 
ues are 6.24 and 5.24, and 6.22 and 5.22 respectively. 

The tidal evolution solution that we presented previ- 
ously would seem to support the hypothesis that the star is 
old, with stellar spin up accounting for the rotation period. 
But the results in Table[5] would seem to add weight to the 
idea that the system is actually younger than expected. The 
likelihood of observing one system in such a condition is low, 
but to observe two such systems as Table[5] implies that we 
have done with WASPs -4 and -19, seems incredulous. 

Considering the remaining lifetimes in the context of 
the total planetary lifetime suggests a different scenario. [13] 
can also be used to calculate the total planetary lifetime 
if the initial semi-major axis is known. We calculated the 
total planetary lifetimes assuming that the value adopted for 
WASP-19 b, ao — 0.0317 AU applied unilaterally across our 
sample, and in Fig.|9] plot remaining lifetime as a function 
of total lifetime. It is immediately apparent that WASP- 
19 b is a special case, being clearly separated from the rest 
of the sample. The remaining lifetime of WASP-19 b is only 
between 1 and 2 percent of the total lifetime, whilst for all 
of the other planets in the sample this figure is greater than 

10 percent. The short remaining lifetime of WASP-4 b is thus 
somewhat misleading, as it represents a significant portion 
(« 12.6 percent) of the total lifetime of the planet. 

There are several ways in which we can analyse the re- 
maining lifetimes calculated here, one of which is to examine 
the probability that we have managed to observe the system 
in its present configuration. We calculated this for the sys- 
tems in Table[5j disregarding those with no literature age 
estimate. For WASP-19 we calculate the probability using 
both our own age adopted age estimate, and the younger age 
of 0.6 ± 0.5 that can be found in the literature. Fig. [lO] dis- 
plays the results for a range of \og{Q'^) values consistent with 




ifetime /Gyr 

Figure 9. Remaining lifetime as a function of total planetary 
lifetime for the sample of transiting planets in Table[5] as cal- 
culated using |13| assuming that the values of log(Q^) and ao 
adopted for the WASP-19 system from our MCMC results ap- 
ply universally. We plot the results for log^Q'^) = 5.52 (squares), 
log(Qs) = 6.47 (circles) and \og{Q'^) = 7.14 (triangles). The data 
for WASP-19 b are denoted by filled symbols. The lines represent 
tremain = 1 percent (red, dash), 5 percent (green, dot-dash) and 
10 percent (dark blue, dot) of the total planetary lifetime. WASP- 
19 b is clearly separate from the rest of the sample, indicating that 
it is in a unique situation. 



the 1 — (T limits derived from the MCMC posterior probabil- 
ity distributions. For the value of log{Q's) adopted from the 
MCMC solution we find that two systems, including WASP- 
19 at our older age estimate, have a probability of observa- 
tion of less than 1 percent. Using the younger age estimate 
for WASP-19 pushes the probability up to 2 percent, which is 



more plausible (see analysis in Hellier et al. (2011 1). Increas- 



ing the value of log(Q!,) to our upper limit of 7.14 increases 
the probability of observing WASP-19 b to approximately 
3.5 percent for the older age, and approximately 9 percent 
for the younger age. However if we consider the lower limit 
of log(Qs) — 5.52 then four systems, including both ages of 
WASP-19, have observation probabilities less than 1 percent. 
This is an entirely implausible situation. 

However this analysis assumes that all of the planets in 
our sample started at the same distance from their respec- 
tive host stars, and that they experience tidal interactions 
of the same strength. It seems somewhat unlikely that this 
accurately represents reality, as the planetary systems in our 
sample show significant variation in their properties. Previ- 
ous studies ( Matsumura, Peale fc Rasio|2010" Hansen|2010 \ 
have found that different planetary systems are likely to ex- 
perience different strengths of stellar and planetary tide, 
so describing these systems with a single value is almost 
certainly unphysical. This would, of course, mean that the 
observational probabilities for several of these hot Jupiters 
could be significantly greater. 

To further muddy the waters we add evidence for an old 
age by carrying out additional stellar model fits. We used the 



updated Padova models (Marigo et al. 20081, which imply 



an age of 7.94^2 93 Gyr, and the latest version of the Yonsei- 
Yale isochrones ( [Demarque et al.|2004[ ) , which return an age 
of 7^3 Gyr. Both estimates agree with the earlier isochrone 
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Table 5. A comparison of the WASP-19 system to a sample of transiting hot Jupiter systems with Mp > Mj in close orbits 
around stars similar to, or cooler than, WASP-19 A. We calculate the remaining lifetime for each planet using equation |13| 
assuming that the stellar tidal quality factor of log(Q'j,) = 6.47 adopted for WASP-19 A is applicable across the entire sample. 
Note that the remaining lifetime for WASP-19 b quoted here does not agree with the value derived from the evolutionary 
track displayed in Fig. [7] 



System Mp /Mj Rp /Rj /M© Rs /R© a /AU P /days tremain /Gyr Age /Gyr 



WASP-19 


1.14 


1.28 


0.95 


0.93 


0.0164 


0.7888399 


0.012 


1 60+2-*-* 
l.DU_(j 79 


WASP-2 


0.847 


1.079 


0.84 


0.834 


0.03138 


2.1522254 


1.708 


1 1 q+*-l a 


WASP-4 


1.12 


1.416 


0.93 


1.365 


0.023 


1.3382282 


0.017 


7 f\~^'^-^ a 
'■"-i.S 


WASP-5 


1.63 


1.171 


1.00 


1.15 


0.02729 


1.6284246 


0.080 


3.0 ± 1.4 ^ 


WASP-10 


3.06 


1.08 


0.71 


0.783 


0.0371 


3.0927616 


1.818 


0.8 ±0.2 " 


CoRoT-1 


1.03 


1.49 


0.95 


1.11 


0.0254 


1.5089557 


0.093 




CoRoT-2 
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0.0281 
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0.0229 


1.4324772 
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Figure 10. The probability of observing the systems in Table[5] 
in their present configuration, as a function of log(Q'g) within the 
range encompassed by the 1 — a errors from the MCMC results. 
Systems without a stellar age in the literature were disregarded. 
The data are offset slightly from the true value of log(Q^) to 
aid clarity for the error bars. Systems for which P > 1 are plot- 
ted at P = 1. WASP-19 is denoted by filled symbols; triangles 
for the probability using the age estimate from the MCMC re- 
sults, squares for the probability calculated using the estimate 
of 0.6 ± 0.5 Gyr prevalent in the literature. At log(Q^) = 6.47 
two systems, including WASP-19 with its older age, show a prob- 
ability of observation that is less than 1 percent. Increasing the 
stellar tidal quality factor increases the probability that we have 
observed the systems in the state implied by their remaining life- 
times. 



fit of |H10[ We also investigated the values of Tef / and ps re- 
quired for a stellar age of 0.6 ± 0.5 Gyr in tfiese two models, 
assuming tiie same metallicity; the results are given in Ta- 
ble|6] The correction required to the effective temperature 
if the stellar density is correct is not improbable, but the 
converse is not necessarily true. Modifying the stellar den- 



sity requires a change to either the stellar mass or radius, 
or both, either of which could have a dramatic effect on our 
estimates of the planetary parameters. This is notwithstand- 
ing the more likely case, which is that both the temperature 
and density are in need of a slight adjustment. A younger 
age for the star does not seem, therefore, outwith the realms 
of possibility. 

For completeness we also attempted to force an age of 
0.6 Gyr by adjusting the stellar metallicity whilst leaving the 
effective temperature and stellar density at the values from 
[HlOl We found that the upper limit oi Z = 0.03 imposed 
on the Padova isochrones prevented us from reaching such 
a low age, leaving us with a lowest age estimate of 3.98 ± 
3.96 Gyr, which could be consistent with a star the age of the 
Hyades. When fitting to the Yonsei-Yale isochrones we found 
that a value of [Fe/H] — 0.32 gave an age of 2.5l'J g Gyr, 
just consistent with the younger age that we were looking 
for. Fitting to the 0.6 Gyr isochrone itself would require a 
much greater metallicity, but given that this value of the iron 
abundance is already 0.05 dex greater than the maximum 
elemental abundance found by HIO (specifically the upper 
1 — a limit on the Calcium abundance), it is clear that such 
a fit would be utterly inconsistent with the current spectral 
analysis. 

It therefore seems that a combination of slightly in- 
creased elemental iron abundance, greater stellar density 
and higher effective temperature could produce stellar model 
fits more consistent with the young stellar age suggested by 
gyrochronology. Indeed if one is using the Padova models 
then such a combination lies just within the existing 1 — cr 
ranges for the three parameters. 

As noted previously, the constraint placed on the age 
by the lithium abundance has the potential to rule out the 
young stellar age. The detection of a substantial lithium 
abundance in a stellar spectrum strongly implies a young 
age for the star, even one that hosts planets, as recent work 
by |Baumann et ar] ( |2010] ) shows that there is no correlation 
between the presence of planets and reduced stellar lithium 
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Table 6. The age estimates obtained through fitting the stellar 
parameters of WASP-19 A to stellar models. Adjusting the stellar 
density, effective temperature or iron abundance in isolation low- 
ers the age that the model fit returns. Obtaining an age of 0.6 Gyr 
required the parameters to be increased beyond their I — a limits 
when done in isolation. If adjusted as a set, the required age could 
be obtained with more plausible values. 



Model 



ps/pe 



/K 



[Fe/H] age /Gyr 



Padova 
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13 ± 





12 


5500 ± 


100 





02 


7 Q4+4-65 
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45 ± 
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5500 ± 
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n (^'^^ +3.35 


Padova 
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13 ± 
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100 
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5500 ± 


100 
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3.98 ± 3.96 
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100 
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5500 ± 


100 
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100 





02 


"■""-0.54 


Yonsei-Yale 
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100 





02 


^•"'-'-0.54 
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5500 ± 


100 
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Yonsei-Yale 
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30 ± 





10 


5625 ± 


100 





12 
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content. The original analysis in |H10| utilised only 34 spec- 
tra from the CORALIE spectrograph ( Queloz et al.|[2000 |, 
but thanks to an ongoing radial velocity observation pro- 
gram we now have access to additional data on WASP-19 A. 



Since the publication of |H10| a further 3 spectra have been 
taken using CORALIE, and 36 spectral measurements have 
been obtained using the HARPS high precision echelle spec- 
trograph (Mayor et al. 20031. We co-added the individual 



spectra into a single spectrum with a higher signal-to-noise 
ratio, from which we were able to improve the constraint on 
the lithium abundance. Hints of a lithium line were present 
but only at the level of the spectral noise, limiting us to an 
upper limit of log(eLi) < 0.5. This allows us to place a lower 
limit on the age of WASP-19 A of 2.0 Gyr (Scstito fc Randich 



2005), making an older age for the star more palatable. 

be considered is the 
Using the ob- 



The final piece of evidence to 
proximity of WASP- 19, to the Roche limit, 
served system parameters and their associated uncertain- 
ties with the formulation of |Eggleton (19831, we calcu- 
late that the Roche limit for WASP-19 b is at 0.0148 ± 
0.0025 AU, in agreement with the observed or bital sepa- 
ration of 0.0164t|3 [;5j^g AU. [HelUer et al] |201l| place the 
planet slightly further from the star at 1.21 times the Roche 
tidal radius, but we note that they have used a different 
formulation for calculating the limit. 

We conclude that the evidence is in favour of WASP- 
19 A being old. Whilst there is evidence for the star being 
young, it is more circumstantial than that which points to an 
older star. The upper limit on the lithium abundance and the 
results from stellar model fitting in particular point towards 
an age in excess of 1.0 Gyr. If WASP-19 A is indeed old, 
then the exploration of tide-governed evolution presented 
herein suggests that WASP-19 b has spun up its host star, 
and might be in the final stages of spiralling into the Roche 
limit. During our exploration of possible evolutionary his- 
tories we found that those which returned a stellar age of 
> 1.0 Gyr tended to exhibit a very short remaining lifetime. 



There were, however, some histories that married an older 
age for the star to a long remaining lifetime for the planet, 
so we are unable to completely rule out that scenario. 



6 SUMMARY AND CONCLUSIONS 

We have calculated an age for the WASP-18 system of 
0.57910.250 Gyr, in agreement wi th th e 0.630lo:53o Gyr age 
found from stellar isochrones by |H09| Using an MCMC al- 
gorithm we find tidal quality factors of log(Qs) = 8.2llo g2 
7.77li:^g. Our results imply that WASP-18 b 
will reach the Roche limit in 0.0761q o44 Gyr, and that in 



and log{Q'p) 



most cases it will cause its host to spin up as it does so. We 
are unable to place stronger constraints the status of the 
system with respect to planetary infall owing to the range 
of evolutionary histories that we find to be compatible with 
its observed parameters, but a large number of the evolu- 
tionary tracks investigated imply that the planet is in the 
process of spiralling in to its host star. Our results for the 
WASP-18 system seem roughly consistent with the work of 



Matsumura, Peale & Rasio (20101 



For the WASP-19 system we found tidal quality factors 

These val- 

+2.8 ' 
0.79 ' 



7+0.67 
'-0.95 

ues give a stellar age of 1.60 



of log(g^) = 6A7t°oil and log(Q;) = 6.75 



Gyr, in broad agreement 
with the constraint of age ^ 1.0 Gyr found by |H10[ and im- 
ply a remaining lifetime of 0.0067^ J Jogj Gyr. We have in- 
vestigated the possibility that WASP-19 A is younger than 
previously estimated, in line with the predictions of gy- 
rochronology. After considering the evidence for both the 
old and young stellar age possibilities, we conclude that 
the older age is more probable based on a reanalysis of the 
spectral lithium abundance and updated stellar model fits. 
We therefore suggest that WASP-19 b could be in the final 
stages of its spiral-in, and could be on the verge of reach- 
ing the Roche limit. We found that this scenario was more 
prevalent amongst those evolutionary histories with a stellar 
age > 1.0 Gyr, but that there were some instances in which 
the older stellar age coincided with a substantial remaining 
planetary lifetime. We are therefore unable to rule out the 
possibility that it will be some time before the planet falls 
into the star. 

Tidal interactions between these two hot Jupiters and 
their host stars will dramatically affect the evolution of the 
stellar rotation periods, counteracting and then reversing 
the spin-down that is expected from evolution according to 
a Skumanich-type magnetic braking law. The observed ro- 
tation periods are irreconcilable with such an evolution of 
the rotation period, and strongly suggest that falling in hot 
Jupiters cause their host stars to spin up during their in- 
ward, tidal interaction governed migration. It is worth not- 
ing that the our results seem to point toward a more di- 
verse range of stellar tidal dissipation strength than is com- 
monly considered in the literature. The ranges of log{Q's) 
that we have attributed to the two systems investigated 
herein are slightly disparate, for which there may be sev- 
eral possible explanations. We turn to the work of |Pinson^ 
[neault, DePoy fc Coffee| ( |2001[ ), which shows that the mass 
of the convective zone is a function of Te//, and therefore 
of spectral type. From their fig. 1, we note that the effective 
temperature of WASP-18 A lies close to the point at which 
the mass of the convective zone becomes negligible, imply- 
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ing a convective zone mass of Mcz ~ 0.001 Msun- WASP- 
19 A, with its much lower effective temperature, would have 
Mcz ~ 0.030 Msun- In our work we have assumed that 
the star rotates as a single body, but we have made no as- 
sumption about where the majority of tidal dissipation takes 
place. If this process occurs in the convective zone then the 
discrepancy in the masses, and therefore depths, of the con- 
vective zones of the two stars could provide an explanation 
for the disagreement as to the value of Q'g, with a larger 
convective zone allowing for more efficient dissipation and 
hence a smaller quality factor. 

This study also provides a warning against using gy- 
rochronology to estimate the ages of hot Jupiter host stars. 
Owing to the tidal spin-up of its host star by the in-falling 
planet, the age that we have found for the WASP-19 plane- 
tary system is greater than the age found using gyrochronol- 
ogy alone. The situation for the WASP-18 system is less clear 
cut, but there is no doubt that gyrochronology will not pro- 
vide an accurate estimate of the system age at all points dur- 
ing its evolution. We therefore suggest that care should be 
taken when applying gyrochronology to hot Jupiter systems. 
The two systems studied herein are extreme examples of this 
class of planet, with extremely short orbital periods and very 
close orbits. It is likely that there will be a critical point in 
semi-major axis space, beyond which tidal interactions will 
be sufficiently weak so as not to influence the stellar rota- 
tion period. For systems in this region gyrochronology may 
well work very well. But for systems such as WASP-18 and 
WASP-19 which have a < acrit, gyrochronology should be 
applied with care. Determination of this critical point will 
require further work, and it is likely to depend on a host of 
factors. 



Lanza ( 2010 1 also suggests that gyrochronology may not 



always provide accurate age estimates for exoplanetary sys- 
tems, finding that plotting Prott~" as a function of Te// 
for planet hosting stars gives a poor fit to the period-colour 
relation of |Barnes| ( |2007[ ). This implies that exoplanet host 
stars are systematically faster rotators than stars with sim- 
ilar ages and properties that do not appear to have any 
associated planets. |Lanza| also finds that the rotation pe- 
riod evolution of F- and G-type planet hosts does not differ 
substantially from similar stars without hot Jupiters, a con- 
clusion with which we disagree although we note that the 
final spiral-in of the planet does not appear to have been 
considered. 

The substantially reduced rotation period that results 
from tidal spin-up may provide a means of detecting stars 
that have either been planet hosts in the past, or that have 
unseen planetary companions that are in the process of spi- 
raling in to the Roche limit. Measurement of such an anoma- 
lous rotation period would provide a strong indication of the 
current or previous existence of a hot Jupiter around that 
star. Searching for such systems could help to pinpoint tar- 
gets for extra-solar planet searches. 
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